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Abstract 

In the presence of ATP, kinesin proceeds along the protofilament of microtubule by alternated binding 
of two motor domains on the tubulin binding sites. Since the processivity of kinesin is much higher 
than other motor proteins, it has been speculated that there exists a mechanism for allosteric regulation 
between the two monomers. Recent experiments suggest that ATP binding to the leading head domain 
in kinesin is regulated by the rearward strain built on the neck-linker. We test this hypothesis by 
explicitly modeling a C^-based kinesin structure whose both motor domains are bound on the tubulin 
binding sites. The equilibrium structures of kinesin on the microtubule show disordered and ordered 
neck-linker configurations for the leading and the trailing head, respectively. The comparison of the 
structures between the two heads shows that several native contacts present at the nucleotide binding 
site in the leading head are less intact than those in the binding site of the rear head. The network 
of native contacts obtained from this comparison provides the internal tension propagation pathway, 
which leads to the disruption of the nucleotide binding site in the leading head. Also, using an argument 
based on polymer theory, we estimate the internal tension built on the neck-linker to be / ~ (12 — 15) 
pN. Both of these conclusions support the experimental hypothesis. 
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I. INTRODUCTION 



Extensive interest has recently been devoted to the understanding of molecular motors, which 
play pivotal roles in cellular processes by performing mechanical work using the energy-driven 
conformational changes. Kinesin, myosin, Fi-ATPase, GroEL, RNA polymerase, and ribosome 
belong to a group of biological machines that undergoes a series of conformational changes 
during the mechanochemical cycle where the molecular conformation is directly coupled to the 
chemical state of the ligand. Although substantial progress has been achieved in understanding 
the underlying physical principles that govern molecular motors during the last decade, major 
issues still remain to be resolved. Specifically, some of the outstanding questions are as follows: 
(a) How is the chemical energy converted into mechanical work? (b) How is the directionality of 
the molecular movement determined? (c) How is the molecular movement coordinated or regu- 
lated? Several biochemical experiments have quantified the kinetic steps , single molecule 
experiments usingoptical tweezers have measured the mechanical response of individual molec- 
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however, are not sufficient to fully address all the questions above. For example, little is known 
not only about the structural details of each chemical state but also about the kinetic pathways 
connecting them. Hence, if feasible, a computational strategy using the coordinates from X-ray 



and/or NMR structures can shed 



ight on the allosteric dynamics of molecular motors. Although 
isl . [igI . \v\ have proceeded towards addressing issue (a) for a 



some initial numerical studies [l4, 
few cases where both open and closed structures are explicitly known, no previous attempt has 
been made to answer issue (c). In this paper we investigate this question in the context of the 
conventional kinesin where the mechanochemical coordination of the motor movement is best 
manifested among the motor proteins. 

One of the experimentally best studied molecular motors is the conventional kinesin (kinesin- 
1) Q, [3], a relatively small sized motor protein that transports cellular materials by walking in 
an unidirectional hand-over-hand manner along the microtubule (MT) filaments. Compared to 
other motor proteins involved in material transport such as myosin and dynein, the conventional 
kinesin has a remarkable processivity, which can travel about a hundred (~8.2 nm) steps with- 
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out being dissociated from the MT. The mechanochemical cycle conjectured from experiments 
suggests that there must be a dynamic coordination between the two motor domains in order 
to achieve such high processivity. The quest to identify the origin of this dynamic coordination 



has drawn extensive attention among the kinesin community. Since Hancock and Howard 20 1 



first hypothesized that the "internal strain" was needed for processivity, the strain-dependent 
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23|. With the aid of 



mechanochemistry became a popular subject in kinesin studies 
optical tweezers, Guydosh and Block recently revisited this issue by monitoring the real-time 
kinesin dynamics in the presence of ATP and ADP • BeFx, a tight binding ATP analog [24]. 
They discovered that, when ADP ■ BeFx was bound to the kinesin, the pause-time of the step 
increased substantially and that the normal step was restored only after the obligatory back- 
step. This suggests that ADP ■ BeFx is released only when the head bound with ADP ■ BeFx 
becomes the leading head (L). Supported by this observation, they advocated a kinetic model 
in which the rearward strain via the neck-linker facilitates the release of the ligand from the 
L 2j]. Stated differently, the binding of the ligand to the L is inhibited because the rearward 
strain constitutes an unfavorable environment for the ATP binding sites of the L. In the present 
study, we focus on the elucidation of the structural origin of the coordinated motion in kinesin 
by adopting a simple computational strategy. 

Better straightforward evidence of the regulation on the nucleotide binding site can be 
obtained when a structure in which both kinesin heads are simultaneously bound to the MT 
binding site is determined. Such a structure will allow us to identify the structural differences 
between the leading (L) and the trailing (T) head. To date, this structure, however, has not 
yet been reported. The only available structures include an isolated kinesin- 1 without the MT 



61], an isolated single-headed kinesin-like (KIFIA ) w ith various ligand states [T\, and a single 



KIFIA bound to the tubulin-dimer binding site [25|]. Therefore, we utilized existing Protein 
Data Bank (PDB) structures and manually built a model system of the two-headed kinesin 
molecule with both heads bound to the tubulin binding sites (see Fig. 2 and legend). This model 
was used to generate an ensemble of structures via simulations. A direct comparison between 
the L and T equilibrium structures shows that the tension built on the neck-linker induces 
the disruption of the nucleotide binding site of the L, which directly supports inferences from 



experimental observations 
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II. RESULTS AND DISCUSSION 



Mechanochemical cycle of kinesin: We begin by reviewing the mechanochemical cycle 
of kinesin molecule on the MT to clarify the importance of dynamic coordination between the 
two motor domains for kinesin processivity. Recent experiments using laser optical tweezers 
(LOT), cryo-electron microscopy, electron pararnagnetic resonance, and FRET, as well as the 



crystal structures at various states 
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28| provide glimpses into the structural and 



dynamical details of how the kinesin molecule walks on the microbubule filaments. Depending 
on the nucleotide state at the binding site, both the motor domain structure and the binding 
interface between kinesin and MT are affected. In particular, a minor change of the motor 
domain coupled to the nucleotide is amplified to a substantial conformational change of the neck 
linker between the ordered and the disordered state. Experimental studies strongly suggest the 



mechanochemical cycle shown in Fig.l 



28l | . The mechanical stepping cycle of kinesin initiates 



with the binding of ATP to the empty kinesin head strongly bound to the MT [{i) ^ (^0]- 
Docking of the neck linker to the neck-linker binding motif on the leading head (X in Fig.l) 
propels the trailing head (Y in Fig.l) in the (+)-direction of MT, which leads to an 8 nm 
mechanical step [{ii) —>■ {ii')]- The interaction with the MT facilitates the dissociation of ADP 
from the catalytic core of L [{ii') — > (iii)]. ATP is hydrolyzed and produces ADP ■ Pi state for 
the T [{iii) — ^ (iv)]. When Pi is released and the trailing head is unbound from the MT, the 
half-cycle is completed [{iv) {{)]. The mechanical step is achieved in a hand-over-hand fashion 



291, 3^. High 



by alternating the binding of the two motor domains (X and Y in Fig.l) to the MT 
processivity of the kinesin requires this kinetic cycle to be stable (remain within the yellow box 
in Fig.l). A premature binding of the ATP to the leading head in the state of {E : MT) should 
be prevented, i.e., the condition fcjf )[ATP]/(fc^™^ + 42) and fcjf ^ [ATP]/(fc^™^ + 42) ^ 
should be satisfied in Fig.l (see Supporting Information for the master equation describing 
the kinetic cycle). ATP binding to the (iii) or (iv) states can destroy the mechanochemical 
cycle of the kinesin. The binding of ATP on the leading head should be suppressed before the 
7-Pj is released from the T. Otherwise, both heads become ADP-bound states, which have a 
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weak binding affinity to the MT, and that leads to dissociation from the MT. Since the kinesin 
has a high processivity compared to other molecular motors, effective communication is re- 
quired between the two heads regarding the chemical state of each of the partner motor domains. 

Two-headed kinesin bound to the microtubule: In the absence of interactions with 
the MT, the individual kinesin monomers fold into identical conformations. To achieve its 
biological function, however, folding into the native structure alone is not sufficient. Coupled 
with the nucleotide and the MT, the two kinesin monomers in the dimeric complex need 
to alternate the acquisition of the native structure in a time-coordinated fashion for the 
uni-directional movement. The currently available three-dimensional structure (PDB id: 3kin, 
structure 2 in Fig.2-A), in which each monomer is in its native state, does not provide such a 
dynamic picture by failing in fulfilling the geometrical requirement of simultaneous bindings of 
both motor domain to the adjacent tubulin binding sites that have an 8-nm gap. The inspection 
of 3-D structure suggests that a substantial increase of the distance between the two motor 
domains can be gained by breaking a few contacts associated with the neck-linker (/39, /310) and 
the neck-linker binding site on the motor domain {(37). To this end, we manipulated the 3-D 
structure of 3kin around the neck-linker of the L and created a temporary structure whose two 
heads bind to the MT binding sites simultaneously. Both L and T have energetic biases towards 
the identical native fold but the interactions with the tubulin binding sites adapt the dimeric 
kinesin structure into a different minimum structure, which is not known a priori. We performed 
simulations (see Methods) to relax this initial structure and to establish the thermal equihbrium 
ensemble of the kinesin molecule on MT (sec Fig.3-A). Transient dimeric kinesin conformations 
corresponding to the steps (iii) and (iv) during the cycle (Fig.l) allow us to investigate the 
structural deviation between L and T of kinesin molecule. This simple computational exercise 
can confirm or dismiss the experimental conjecture regarding whether the mechanochemical 
strain significantly induces regulation on the nucleotide binding site and also if it occurs in the L. 

Catalytic core of the leading head is less native-like on the MT: Since the nucleotide 
binding and release dynamics is sensitively controlled by the kinesin structure, we assume that 
the nucleotide molecule has an optimal binding affinity to the kinesin motor domain in the 
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native structure. For function there is a need to understand how the native structure of the 
kinesin motor domain is perturbed under the different topological constraints imposed on the 
dimeric kinesin configuration by interacting with the MT. The equilibrium ensemble of the 
structures shows that the neck-linker is in the docked state for the T but undocked for L. In 
comparison to the native structure, the overall shape of the nucleotide binding pocket in the T 
is more preserved. As long as the MT constrains the two heads 8-nm apart, this configuration 
is dominant in the thermal ensemble (see Fig.3-A). 

Global shape comparison : There are in principle multiple ways to quantitatively compare the 
two motor domain structures. To assess the structural differences, the radius of gyration {Rg — 
l/2N'^'^-j(Ri — RjY) of the two motor domain structures from the equilibrium ensemble are 
computed (see Fig.3-B). Because the neck-hnker and the neck-helix adopt different configurations 
relative to the motor domain in each monomer, we perform a Rg analysis for the motor domains 
only (residue 2-324). The Rg distributions show that the L is slightly bigger than the T both in 
the size {{{Rg){L) - {Rg){T)) ~ 0.4 A) and in the dispersion {{ul - (Tt) ~ 0.05 A). Meanwhile 
the Rg for the native state (3kin) is i?^"**^^ — 19.4 A. Clearly, the sizes of both of the heads in 
the thermal ensemble are expanded at T = 300ir as compared to the native structure. The size 
alone does not tell much about the difference between the structures. 

The RMSD relative to the native structure and between the two motor domains 
(residues 2-324) computed over the equihbrium ensemble gives Y{MSD{T\native) — 2.0 A, 
RMSD(L|natiT;e) = 9.4 A, RMSD(r|L) = 9.6 A, where RMSD(A:|y) is the RMSD between 
conformations X and Y. If the q;6 helix is excluded from the RMSD calculation of motor 
domain (residues 2-315), then RMSD(T|native) ^ 1.8 A, RMSD(L|native) ^ 3.8 A, and 
RMSD(T|L) 3.9 A. The RMSD analysis shows that the q;6 helix significantly contributes 
more for the deviation of the leading head from its native state than the T. 

Additional detailed comparisons with respect to the native state can be made using the 
structural overlap function oi i — j pair, Xiji which is defined as Xij = {^{Rij — Rlj)) where 
5{Rij — R°j) = 1 if \Rij — R°j\ < Ttoh otherwise 5{Rij — R°j) — 0. R°j is the distance oi i — j pair 
in native state, where rtoi = 1 A. By setting R°j values identical in both heads (i.e., both heads 
have the same native state), we compute the Xij values for the trailing and the leading heads. 



6 



respectively. The relative difference of the Xij value between T and L, Xij is defined by 



X^J{T)-X^J{L) 



[ {X^AT)=0) 

which quantitatively measures the structural difference of the two heads. Based on the Xij value 
(Fig.3-C), the distances between the MT binding motif of the T (Lll, L12, a4, a5) and other 
secondary structure units aO, (32, al, a2, a6) are 50 % more native-like than in the L. 

Conserved native contacts in trailing head reveals the strain propagation pathway in leading 
head: A direct measure of similarity to the native structure is the fraction of native contacts 



preserved in the thermal ensemble 31|- Since we assume that ATP affinity is optimized in the 
native state, we can readily assess the quality of the structure using this measure. We quantify 
the nativeness of a pair using qij{^) = (0(-Rc — Rij)^ij)i where Ajj = 1 if j residues are in 
contact at the native state < = 8 A), and Ajj = otherwise. (with ^ = T or L) 

is obtained by averaging over the thermal ensemble. When qij{C) is averaged over all the native 
pairs, the average fraction of native contacts, {Q) is calculated as {Q){C) = ^/^Q^i<j%iO 
where Nq is the total number of native pairs. For the T and L conformations, {Q){T) = 0.86 
and {Q){L) = 0.82, respectively. The relative difference of native contacts between the two 
kinesin heads at the pair level, Qij, is quantified similarly to EqUJas 

f (q (T) ^ 0) 

\ (g.,(T) = 0) 

In Fig. 4, Qij is color-coded based on its value. As expected from the equilibrium ensemble, 
conspicuous differences are found around the structural motifs having direct contacts with 
neck-linker, giving Qij > 0.5. Quantitative inspection of the other contacts is illustrated in the 
structure. We color the kinesin head structure based on the Qij value. The residue pairs are 
colored in magenta if 0.2 < Qij < 0.5, red if Qij > 0.5, where the positive Qij signifies that the 
native contacts in trailing head are more intact. The residue pairs are colored in light-blue if 
—0.5 < Qij < —0.2, blue if Qij < —0.5. More intact contacts, when the trailing and the leading 
head are compared, are visualized by yellow line in Fig.5-B. Our analysis not only shows that 
there is higher probability of the formation of native contacts present in the T in comparison 
to the L, but also suggests how the tension is propagated towards the nucleotide binding site 



to disrupt the nativeness of the nucleotide binding pocket in the leading head. As expected, a 
dense network of intact contacts are found between the neck-linker [jSlO) and the neck-linker 
binding motif {(37). This network continues along the a6 helix, perturbing L2, (31, aA, and 
finally reaches the nucleotide binding site (see SI Fig. 6 for the nomencaltures of the secondary 
structures). It is surprising that the disruptions of native contacts are found particularly 
in the nucleotide binding site, which is believed to be the trigger point for the allosteric 
transition. All the important nucleotide binding motifs (P-loop, switch-1, switch-2, and N4) are 
recognized by our simulational analysis using a nonlinear-Hamiltonian (see Supporting Informa- 
tion for the comparison with linear-harmonic potential represented as Gaussian network model). 



Estimate of the tension in the neck-linker: The deformation of the leading motor 
domain is caused by the internal tension in the neck linker. The tension on the neck linker is 
estimated using the force (/) versus extension (x) relationship of a worm-like chain model [32], 

1 X 1" 
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bp 



(3) 
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where Ip is the persistence length of the polymer and L is the contour length. L ^ 5.7nm for 
the 15-amino acid neck-linker (residue from 324 to 338) (= 15 x 0.38nm), and in the equilibrium 
ensemble of structures, x ~ 3.1 ± 0.8 nm. Assuming that Ip ~ (0.4 — 0.5) nm [33] for this 
segment, we estimate a tension / = 12 — 15 pN. By integrating EqJS] for (0 — 3.1) nm. The 
tensional energy stored in the neck-linker is obtained, which is 17 pN ■ nm ^ AksT. About 20% 
of the ATP hydrolysis energy (~ 25/^^^) is stored in the neck-linker and directly perturbs the 
nucleotide binding site of the L, whereas mechanical action to the T is dissipated through the 
dense network of contacts formed between the neck-linker {(310) and the neck-linker binding site 
{(37). 

For a given extension x{< L), when the length of the neck-linker is varied by 
6L, the variation in the length of the neck-linker can affect the effective tension as 
f{L + 6L) - f{L) ^ -kBT/lp[l/{4{l - xjLf} x (2x/L)/(l - xjL) + x/L\{5L/L). For 
a given 5L/L, the resulting tension change may be significant depending on the value of 
the extension x. Experimentally, the kinesin dynamics has recently been studied by varying 
the linker length by introducing a spacer composed of amino acids [34]. Since 5L = 0.38 
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nm for the insertion of a single amino acid, the lengthening of the linker leads to a re- 
duction of the tension by A/ ^ —2pN. In light of the force values controlling the kinesin 
dynamics in LOT experiments, which is (0 — 7) pN and 7 pN is the stall force, a value 
of A/ ^ —2pN can be significant. According to Hackney et. a/.'s experimental analysis 
the processivity is reduced by ~ 2-fold when a single amino acid is inserted or deleted, and 
6 or 12 additional amino acids resulted in to a 3—4 fold reduction in the kinetic processivity [34] . 



III. CONCLUSIONS 



Because of the size and time scale spanned by a typical molecular motor as well as the 



lack of crystal structures, theoretical study 



such as master equation descriptions 
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3ased on the structure is not a common approach 
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37H or Brownian ratchet models 



Knowledge 



of structural details, however, is the key ingredient to understand the mechanochemistry of 
molecular motors. In the present study, we propose computational strategies to resolve this 
problem in kinesin dynamics, particularly the nucleotide binding regulation mechanism between 
the two motor domains. By building a kinesin model on the tubulin filament, we explicitly 
identified the effect of internal tension (/ = 12 — 15 pN) on the front kinesin head domain 
and showed that the tension propagation to the leading head provokes the switch-related motifs 
(P-loop, switch-1, switch-2) in the nucleotide binding pocket. Assuming "the nativeness as 
a criterion for the optimal nucleotide binding condition", we concluded that the nucleotide 



binding pocket in the leading head is not favorable (or partially unfolded or cracked 391]) for 



the nucleotide binding while the trailing head is bound to the MT. This conclusion ex plai ns the 
recent real-time single molecule traces of kinesin generated by Guydosh and Block 24!]. The 
reduction of ligand affinity of the leading head due to the rearward tension benefits the high 
processivity in two ways. First, the premature ATP binding is inhibited before the chemical 
reaction (ATP hydrolysis. Pi release) in the trailing head is completed, during which the trailing 
head is tightly bound to the MT and the tension is built on the neck linker. Second, the release 
of ADP is facilitated in the leading head, which accelerates step {ii') [iii) in Fig.l. The 
latter point is consistent with Uemura and Ishiwata's experimental observation of a seven-fold 
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increase of the ADP dissociation constant from the monomeric kinesin head in the presence of 
rearward loading 21 1. 



It is noteworthy that global measures characterizing the two head domains do not distinguish 
the qualitative differences between the two heads. Major changes in the strained region such 
as the kinesin neck-linker may be associated to minor global differences in the motor domain. 
Exactly how these minor global differences that are localized in a small volume is amplified in 
the process of allosteric transition is the another key issue in understanding molecular motors. 

The internal tension regulates the interaction between the kinesin and the nucleotide. Con- 
versely, the interaction between kinesin, nucleotide, and MT also switches on and off the internal 
tension. For conventional kinesin, mechanical and chemical mechanism are closely correlated, 
producing a remarkable processivity of the kinesin movement on the MT. 

IV. METHODS 

The simulations were performed using two classes of energy function. One is the standard 



structure based (SB) potential [40|] and the other is the self-organized polymer (SOP) potential 
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42| (see Supporting Information and SI Figs. 7-8). For these two different topology- 
based potential functions, we obtained qualitatively identical results. Results from the SB 
potential are presented in Figs. 3-5 while results from the SOP potential are in SI Fig. 9 in 
Supporting materials. This suggests that the results are robust as long as the information of 
native topology with a nonlinear form of energy potential is used as an input. 

Energy function— Structure based potential: Using the structure in Fig.2-C as a 
starting structure, we simulated and sampled an ensemble of two-headed kinesin configurations 



on the MT. We performed the Langevin simulations of the SB model 



401 1 whose equation of the 



motion of each interaction center is integrated by a Verlet algorithm. The energy potential is 
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given as 
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where the energy Hamiltonian is divided into intramolecular interactions for the kinesin 
molecule and intermolecular interactions at the kinesin-MT interface. The superscripts K and 
K — tuh denote the kinesin and the kinesin-tubulin interaction, respectively. Because our focus is 
on the kinesin dynamics, we fixed the coordinates of the MT in space. Since the length scale of 
the kinesin geometry is small (< lOnm) compared to that of MT (diameter ~ 24nm, persistence 
length ~ 1mm), the explicit computation of the dynamics of the entire MT structure, in which 
the 13 protofilaments constitute the cylindrical geometry, should not qualitatively change our 
conclusions. The first and the second term define the backbone interactions. The bond distance 
rj^i+i between the neighboring residues i and i + 1 are harmonically constrained with respect to 
the bond distance in native state r°i^-^ with a strength Kj. — 20 kcal/ {mol ■ A^). In the second 
term, the angle 9 is formed between residues i, i + and i + 2 with Ke — 20 kcal /{mol x rad?). 
9° is the angle of the native state. The third term is the dihedral angle potential with 
K^^^ = l.Okcal/mol and K^^ = 0.5kcal/mol that describes the ease of rotation around the 
angle formed between successive residues from i to i + 3 along the backbone. The Lennard- Jones 
potential is used to account for the interactions that stabilize the native topology. A native 
contact is defined from the pair of interaction centers whose distance is less than = 8 A 
in native state for |i — j| > 2. If i and j residues are in contact in the native state, Aij = 1, 
otherwise Ajj = 0. Aided by Aij we assign stabilizing potential for native pairs and repulsive 
potential for non-native pairs. We assign e/j = 1.8 kcal /mol for the intra and inter neck-helix 
(residue> 338) interactions to secure the coiled-coil association between the neck-helices. For 
other kinesin residue-residue interactions, we set €h — ei — 1.0 kcal /mol regardless of the 
sequence identity. The parameters determining the native topology r"j and Ajj are determined 
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from the crystal structure of human kinesin (PDB id : Skin) and incorporated to the traihng 
kinesin (T) and the coiled-coil whose structure is shown in Fig.2-C. To constitute an identical 
fold condition, we transferred topological information in trailing head (T) to the leading head 
(L) by substituting rj'^-, 91, 0° and A,^- from T to L, i.e., r1^{L) = r°^{T), e°{L) = e°{T), 
4>°{L) = (f)i{T) and Ajj(L) = Aij{T) for all i and j. The kinesin-tubulin interaction parameters 
(e/j, ei, r°i^, and A*^) are similarly defined as kinesin intramolecular interaction parameters 
except for the slightly larger native contact distance (i?^"*"'' = 10 A). The parameters r?^ and 
A*;;, defining the interface topology between the trailing kinesin head (T) and the tubulin are also 
transfered to the interface topology between the leading head L and the next tubulin binding site. 

Simulations : The initial structure, whose two heads are constrained to be oriented on 
the tubulin binding sites, is relaxed under the SB or SOP-Hamiltonian and subsequently the 
equilibrium ensemble of the structures is collected from the low friction Langevin dynamics 
simulations at T = 300i^'. The position of the interaction center is integrated using 

mf = —Cf — ^— + T (5) 
or 

where ( is friction coefficient, — ^ is the conformation force, and T is the random force 
satisfying {T{t) ■ T{t')) = — t') where the integration time (h) is discretized. In low 

friction Langevin dynamics, natural time is given by tl = {ma^ /thf'^. We chose C = 0.05r^^ 
and h = 0.0025ri. Low friction is deliberately chosen for the purpose of effectively sampling 
the conformational space {a,^. Under such conditions the resulting dynamics as a function of 
time step should not be taken parallel to the real time dynamics. To produce an overdamped 
dynamics it is essential to integrate the motion by neglecting the inertial term as well as 
choosing a high friction coefficient that amounts to the water viscosity (~ IcP). 
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VI. FIGURE LEGENDS 



Figure 1: Mechanochemical cycle of kinesin. The subscripts X and Y refer to each of the 
kinesin head, E denotes the empty head, and (:MT) is appended if the head is strongly bound 
to the MT. 

Figure 2: Procedure to construct the two-headed kinesin/MT-protofilament model. A. 
Three structures from protein data bank (PDB) are used. 1. Single headed kinesin (KIFIA) 
bound on tubuhn (PDB id: liaO). 2. Two-headed kinesin (PDB id: 3kin). 3. Two consecutive 
tubulin complexed to the stathmin-like domain (PDB id: Iffx). We overlaped the chain (A), 
(B) (blue in figure) of Skin onto the chain K of liaO and the chain A of Iffx (a-domain) onto 
the chain (B) of liaO, which leads to the structure in B. The structural homology (Cq, backbone 
RMSD=1.6A) between KIFlA and a head of the two-headed kinesin are sufficient that one of 
the kinesin heads fits to the tubuhn binding site. While the sequence difference between KIFlA 
and the conventional kinesin (sequence identity ~ 45%) may affect the strength of interactions 
between kinesin and tubulin, leading to a different binding affinity of conventional kinesin from 
KIFlA, we assume that the binding orientation of the two-headed kinesin is similar to that 
of KIFlA on the tubuhn. After the structure overlap, the chain (C), (D) of Skin is internally 
rotated around a few positions in the neck linker (S24-SS8) until the chains C, D are placed in 
the vicinity of binding site of tubulin that is designed to be identical to the interface between the 
kinesin rear head and the tubulin. We performed the simulation to relax the kinesin structure 
on the MT and obtained the structure shown in C. 

Figure S: The ensemble of structures and structural comparisons between two heads using 
Rg and Xij. A. The thermal ensemble of structures is illustrated using the multiple structures 
obtained during the simulations. Different colors are used to distinguish the motor domain 
(residues 2— S2S) from the neck-linker and the tail part (residues S24— S70). Substantial vari- 
ations of the neck- linker/tail position in ensemble show its flexibility. The nucleotide binding 
sites in the L and T are indicated by the arrows. On the upper-right corner of the panel are 
two crystal structures of kinesin (Skin). One is the view from the top (left) and the other is the 
view towards the nucleotide binding pocket (right). When compared with the top view of Skin 
crystal structure, it is visually clear that the nucleotide binding pocket is more intact in the T. 
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B. Analysis of kinesin motor domains (2-324) using the radius of gyration (Rg). Histograms of 
Rg collected over the ensemble are fit to a Gaussian distribution. For the traihng head (red), 
{Rg)T = 20.3 A, (Jt = 0.13 A. For the leading head (blue), (i?^)^ = 20.7 A, aL = 0.18 A. C. 
Analysis using the structure overlap function. Xij value is color-coded on the right (see the text 
for the definition of X^). 

Figure 4: Comparison of two heads using the fraction of native contact. A. The average 
contact map for kinesin (left). The relative difference of average contact map between trailing 
head and leading head with respect to trailing head is shown on the right. The red arrows mark 
the set of contact pairs whose Qij value is greater than 0.5 (50%). B. The result of Qij in A is 
redrawn using 3-D plot for clarity. 

Figure 5: The structure of kinesin on MT colored based on the protocol discussed in the 
main text. A. The residues colored in red and magenta on the trailing head (structure on the 
left) are those maintaining more contact than the residues in the leading head (structure on the 
right) . Some of the residues in colored region of the trailing head are involved with the nucleotide 
binding site. B. The enlarged view in cartoon representation from neck-linker docking site (left) 
and from ligand binding site (right). Sphere representations in orange are ADP molecule. The 
network of contact pairs are depicted in yellow lines. 
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Supporting Information 



Structural Details of Kinesin. Detailed knowledge of the structure is fundamental 
to understand the working mechanism of a biological nanomachine. A monomer of the 
conventional kinesin is structurally categorized into three parts : The head (or motor domain) 
(residues 2—323), the neck linker (residues 324—338 : (39, /310), and the neck-helix (residue 
339— : a7). The head region is composed of eight /5-strands flanked with three a-helices [al, 
a2, aS on one side (Fig. 6A), and a4, a5, a6 on the other side (Fig. 6B)] on each side of the 
/9-sheet. One side contains the binding motif (a4, a5, a6, L8, Lll, L12) for the microtubule 
(Fig. 6A) and the other side provides a nucleotide binding pocket (Fig. 6D). The nucleotide 
binding site in the kinesin head region is structurally homologous across the motor protein 
and the G-protein superfamilies. The structural motifs around the nucleotide binding site, 
such as the P-loop (Nl) (86-93), switch-1 (N2) (199-204), switch-2 (N3) (232-237), N4 



14-17) are accordingly designated 4J, |45|. The crystal structures with the different nucleotide 



states suggest that the presence or the absence of 7-Pj is sensed by these motifs and that the 
structural changes of the motifs are related to the allosteric transitions. The ordered state of 
the neck-linker, which is extended from the N terminus of kinesin head and is composed of 
two beta strands (/59 and /310), forms contacts with the N-terminal region of the /57-strand. 
Further extension from the neck-linker leads to the neck- helix (a7- helix), through which two 
monomers form a dimeric complex. 



Computations of Residue Displacement Cross-Correlation Using Elastic Gaussian 
Network Model and Simulation under SB Potential. For given coordinates of a complex 
three-dimensional structure, the dynamical property of an object can be extracted at the zeroth 
order by investigating its topology. The Gaussian network model (GNM) is the simplest possible 



46 



47 



48|. The 



method to study large biomolecules using the corresponding minimal topology 
GNM views the biomolecular construct as a collection of beads connected with harmonic springs 
with strength 7. The connectivity between the beads is purely determined by the cut-off distance 
parameter Rc- For a pair of beads i and j, whose distance Rij{= \Ri — Rj\) satisfles Rij < Rc, 
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the harmonic potential constrains the position of beads via 

^ = E i(4 - Rij^^iRc - Rij) = |<^R^ • r • 5R, (6) 

where ©(■■■) is the Heaviside function, (5R^ = {5Ri,..., SRn) with 5Ri — Ri — R°, and Fj^ = 
\ dSR dSR (^'i) element of Kirchhoff matrix F. Use of the harmonic potential amounts to 

the expansion of the potential H{{R}) at the potential minimum {R°} as iJ({R}) = iJ({R°}) + 

ISK^q£§^6R +■■■. Since the partition function of GNM is given by ^ J L>[(5R]e-^^ = 

/ \ 1-3/2 

det I ^ I , the correlation between the spatial fluctuation of two residues is expressed 



\2TrkBT 

using the inverse of Kirchhoff matrix, 

ZAP AR\ 2kBTd log Zn ^^/p-n 

For i — j, the mean square displacement of the i*^ residue, (SRf), corresponds to the B- factor 
(Debye- Waller temperature factor) as Bi — ^{SR"^). A comparison between the B- factor 
and mean square displacement (MSD) from the GNM determines the effective strength of the 
harmonic potential that stabilizes the structure. Note that the quality of the MSD in GNM is 
solely controlled by the Rc value, thus we scaled the MSD with SksT/^ for GNM analysis. 

We applied the GNM analysis with Rc = 8 A on the two-headed kinesin whose both heads fit 
to the adjacent tubulin binding site, and then computed the cross-correlation matrix as shown 
in Fig.7A. The cross-correlation value Cij, {SRi ■ SRj) scaled by 3kBT/'j, shows that except for 
the neck-helix region the amplitude of correlation in leading head is always larger than that of 
the traihng head. This is expected since the neck-hnker of the leading kinesin is detached from 
the motor domain. The residues in the network with less coordination number are subject to a 
larger fluctuation. The relative difference of the cross-correlation between the leading and the 
trailing kinesin using 

c _ Cij{i,j eL)- Cij{i,j e T) 

C,,(z,jGL) 

is illustrated in Fig. 7B. Fig. 7C shows the auto-correlations (or mean square displacement), 
which are the diagonal elements of the C^j matrix. 

GNM analysis is useful in analyzing the fluctuation dynamics of the stable structure at the 
residue level in the basin of attraction where the basin is modeled as a quadratic potential. 
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However, the expansion of the potential minima up to the quadratic term is justified only if 
the fluctuation 6K is small. The amplitude of fluctuation in biological systems at physiological 
temperatures (T ~ 310K) is most likely to exceed the limit beyond which nonlinear response 
is no longer negligible. In order to take this effect into account, the Hamiltonian should be 
expanded beyond the linear response regime. This procedure indeed reverses the simple idea 
that Tirion and Bahar et al. \^}] have proposed in the context of GNM analysis. However, 
minimal inclusion of the nonlinear term can be useful by increasing the susceptibility of the 
structure. Once the nonlinear term is included, a simple analytical expression such as Eq|7] is 
not available. Thus, we resort to the simulations. 

The analytically obtained quantities, {SRi ■ 6Rj), Sij in Fig. 7 can also be calculated over 
the thermal ensemble of structures obtained from simulations using a nonlinear-Hamiltonian 
(see Fig. 8). The flrst conclusion drawn from the simulational analysis is similar to the GNM 
in that the leading head experiences larger fluctuations. Secondly, the position and relative 
amplitude of the MSD peaks, reproduced using the simulation results, shows a good agreement 
with GNM results. However, the direct comparison of Cij (or 6ij) between Fig. 7 and Fig. 
8 shows that the simulation results from the nonlinear-Hamiltonian display a more sensitive 
pattern of cross-correlations. The pronounced amplitude of Cij (or 6ij) suggests a strong spatial 
correlation between residues i and j. 



Alternative Energy Function : SOP Potential. An alternative potential function for 
the SB potential used in the main text is the self-organized polymer (SOP) potential that was 
recently adopted for simulations of the mechanical unfolding of large molecules of RNA and 



proteins 



41 



42| as well as the allosteric dynamics of GroEL [IG]. The energy Hamiltonian is 
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defined as 



Nk-1 




A,- 



a 



'ik 



(9) 

i = l k = l L\\*'<^/ \ IK / J \ IK / 

The first term is for tfie cliain connectivity of tfie kinesin molecule. Tfie finite extensible 
nonlinear elastic (FENE) potential 4^ is used with k = 2Qkcal / {mol-K^)^ Rq = 2 A, and 
Tj^j+i is the distance between neighboring interaction centers i and i + 1. The Lennard- Jones 
potential interactions stabilize the native topology. A native contact is defined as the pair of 
interaction centers whose distance is less than R^, = 8 A in native state for \i — j\ > 2. If i and j 
sites are in contact in the native state, Ajj = 1, otherwise Aj^ = 0. We used eh = 1.8 kcal/mol 
in the native pairs, and ei = 1 kcal/mol for non-native pairs. To ensure the non-crossing of 
the chain, we used a 6^^ power potential in the repulsion terms and set a = 3.8A, which is 
typical Ca — Ca distance. The parameters determining the native topology, r°- and Ajj, are 
adopted from the trailing kinesin (X) whose structure is shown in Fig. 2C. We transferred the 
topological information in the trailing head (T) to the leading head (L) by substituting r°- 
and Aij from the T to L. Kinesin-tubulin interaction energies are similarly defined as kinesin 
intramolecular interaction energies with slightly different native contact distances. We set the 



cut-off distance for the native interactions between the kinesin and the tubulin as R, 



K-tub 
C 
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A. The parameters, and A*;;., defining the interface topology between the kinesin head T and 
the tubulin is transfered to the kinesin head L and the next tubulin binding site. Using the SOP 
potential, we obtained qualitatively identical results as those obtained from the SB potential. 
The nucleotide binding pocket of the front head is disrupted in the dimeric kinesin configuration 
whose both heads are bound to the tubulin binding sites. The figures corresponding to Fig. 3^ 
and C and Fig. 5 are regenerated using SOP model in Fig. 9. 
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Master Equations for the Mechanochemical Cycle of Kinesin Described in Fig. 

1. In the limit when the dissociation of dimeric kinesin from the microtubule is negligible, the 
kinetic equation describing the dynamic cycle shown in Fig. 1 is written as 



dt 



dP, 



{ii 



dt 
dP(ii' 



dt 

dP(^iii 



dt 



dPi 



(iv 



dt 

dP{iii' 



dt 
dP(iv' 



dt 

dP{v 



dt 



-kbi[ATP]P(^i) + krP(ii) + kdMTP{iv) + kaP{v) 

-{kr + kD)Piii) + kbi[ATP]P^i) 

-kdADpP{u') + kDP(ii) 

-{kh + ki;;'\ATP])P^ui) + kdADpPiii') + kf'^P^iw^ 
-{kdMT + k'i;;\ATP])P^^,,) + khP^^,u) + k^r^PiiV) 

-{k^r^ +k^l)Piiv') + kt:\ATP]P^,^^ 



(10) 



where P(a) is the probability of finding the molecule in a mechanochemical state a{— i, ii, ■ • • v) 
with ^Q-P(a) = 1- The steady state solutions by setting = leads to 
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P(iii) 



kuiATP] 



dt 

1 X 



1 X 



(y) 



Z P -1 



ko Z 



f'dADP 



Z 



Pi- — 



= kdMT 
kdMT 



y 



kh 



_j_ diss 

kh 

1 + pis^j^ijv) 



1 + 



'^diss 
''^dMT 



dMT 



1 + 



+ 



'^diss 
ka 



1 + 



y 



1 + 



diss 
ka 



ku[ATP] V A; 



— j + k^ + k^^^p 



X 



j^m^^ kr{ATP\ ^(,.) ^ fci£[ArP] 



(11) 



diss 
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When the average velocity at steady state is computed using v = d{ki,i[ATP]P(^i) 
X j one can write the velocity in the form of Michaelis-Menten equation. 

fc* \ATP\ 



V = d- 



il+Q{[ATP])) 1 



d- 



u1[ATP] 



1+Q{IATP\) ^ 



+ [ATP] 



"■0 



+ [ATP] Km + [ATP] 



(12) 



where {k*) ^ = k^^ + k^X^p + k^l,^ + k^^ ^ 

1 - 



Q{[ATP]) 



k* 
h 



J (Hi) \ / ■ ■ ■\ 

1 + %^ ] 



1 _l diss T/' 



iiii) 



- 1 



+ 



k 



(IMT 



1 + I 1 + ^ ] i^-) 



, (Hi) / ■ ■ ■\ 



(13) 



and d = 8.2 nm (the gap between neighboring tubulin binding sites). If the dissociation from 
the microtubule is suppressed by small Km^^ and Km^ then Q ^ 0. Depending on the rate 
constant, Q can be either positive or negative. Although the large dissociation constant reduces 
the processivity, the presence of dissociation can increase the effective velocity of kinesin if Q < 0. 
The second and the third expressions following the equality sign are given to compare our result 
with the (N=2)-model of Fisher et. al. 36| and Michaelis-Menten kinetics, respectively. If 
Q = the (N=2)-model analysis on the experimental data by Block and coworkers 3| predicts 
k* = u1 = 108s-\ kbi = k^ = 1.80fxM-h-\ and k*kr/kD = w° = 6.0s-^ (kn = 18 x kr). This 
sets the lower bound for the parameters as ki), kdADP, kdMT, i^h 
The average run length of kinesin, L, is calculated using 



ku > 108 s-\ 



L = dx {I) = dx^lil- P(v)yP( 



d X 



P 



(v) 



p, 



(14) 



^=1 - (-) 

where / is the number of mechanical steps of kinesin. If the probability of the dissociated 
kinesin is small (P(t,) ^ 0), (i.e., if [{ADP)y — {ADP)x](v) is negligible in Fig. 1), then 
L ^ d/P(^yj = v/kdiss where kdiss is the dissociation rate. 
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VII. SI FIGURE CAPTIONS 



Fig. 6: Kinesin structure. (A) Top view wlien kinesin is bound to tubulin. The helices 
on the top with respect to /3-sheet are colored in pink. (B) Bottom view. The helices on the 
bottom with respect to /?-sheet or on the side of the tubulin binding interface are colored 
in lightblue. Lll loop, which is not observable in the crystal structure because of disorder, 
is tentatively drawn with a dashed line. (C) Side view. Note that the neck-linker (/59, /310) 
is connecting the a7 neck-helix with the a6 helix in motor domain. (D) View around the 
nucleotide binding site. P-loop, switch-1, switch-2, and N4 regions, which are relevant to 
nucleotide binding, are colored in red with annotation. 

Fig. 7 : Analysis of the kinesin equilibrium dynamics using the Gaussian network model 
(GNM). {A) The cross-correlation map of the residue fluctuation for the two-headed kinesin 
structure is shown in Fig.2C. (B) The amplitudes are color-coded based on its value. The 
relative difference between the trailing and the leading head with respect to leading head {6ij) 
is plotted on the right. ((7) Mean square displacement of residues in the trailing head (red) and 
the leading head (blue) are plotted on the same plot. A comparison of the amphtudes shows 
that the leading head fluctuates more than the trailing head. The relative difference between 
the two plots with respect to leading head (Su) is plotted on the right panel. 

Fig. 8 : Analysis of the kinesin equilibrium dynamics using an equilibrium ensemble gener- 
ated from the simulations under SB-Hamiltonian. The legends for A- C are identical to SI Fig. 7. 

Fig. 9 : Results of the strain induced regulation in the kinesin dimer on the MT generated 
using the SOP potential (Eq|9]). Comparisons between the results obtained with the SOP 
potential and with the SB potential confirm qualitatively the identical conclusions. 
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FIG. 9: 
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